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ABSTRACT 

Recent observations have reported that some gas-rich dwarf irregular (dirr) 
galaxies appear to have spherical distributions in the outer underlying old and 
intermediate-age stellar populations (e.g., NGC 6822). These observations imply 
that some dirr's have two distinct (or "two-component") structures, i.e., inner 
disky and outer spherical ones, though the number fraction of dirr's with such 
structures remains observationally unclear. We discuss how such two distinct 
structures are formed during dIrr formation based on observations and simula- 
tions. Our numerical simulations show that the remnants of mergers between two 
gas-rich dirr's with initially extended gas disks can have both extended spheroids 
composed of older stellar populations and disks composed mostly of gas and 
young stars. The simulated remnants with two distinct structures can be still 
identified as dirr's owing to the presence of star-forming regions. The structural 
properties of outer spherical structures in dirr's formed from dirr-dirr merging 
depend on initial conditions of merging, which suggests that outer structures in 
dirr's can be diverse. We also discuss other possible physical mechanisms for the 
formation of outer spherical structures composed of older stars in dirr's. 

Subject headings: galaxies: irregular - galaxies: ISM - galaxies: dwarf - galax- 
ies:evolution - galaxies: stellar content 



1. Introduction 



Luminous galaxies like the Milky Way can be formed from merging and accretion of 
numerous low-mass galactic building blocks like dwarf galaxies in hierarchical galaxy forma- 
tion scenarios based on the cold dark matter cosmology (e.g.. White & Rees 1978). Physical 
properties of dwarf galaxies thus have been discussed extensively in many different contexts 
of galaxy formation and evolution, such as physical mechanism for galaxy-scale star forma- 
tion process in dirr's (e.g.. Hunter 1997) and evolutionary links between different types of 
dwarfs (e.g., Grebel et al. 2003). Stellar populations, structures and kinematics of HI gas. 
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and spatial distributions of young stars in gas-rich dirr's particularly have attracted much 
attention from theoretical works on formation of stars and star clusters (e.g., Gallagher & 
Hunter 1984; Hunter 1997) and inner structures of their dark matter halos (e.g., Burkert 
1995). 

Although previous observations revealed fundamental physical properties of star- forming 
regions, young stellar populations, and HI gas in dirr's (e.g.. Hunter 1997), structural and 
kinematical properties of underlying older stellar populations remain fairly unclear. Obser- 
vational studies of dirr's based on broadband imaging {UBV JHK) revealed that the struc- 
tures of main disk components in dirr's can be described as exponential profiles (Hunter 
& Elmegreen 2006). Furthermore observational studies on mid-infrared properties of dirr's 
suggested that the overall structure of a dirr is the same in different passbands for most 
dirr's (Hunter et al. 2006). Since these observations are based on integrated stellar hghts, 
they did not allow the authors to investigate separately structures of stellar populations with 
different ages in detail. 

Recent observational studies on spatial distributions of individual older (e.g. RGB stars) 
and intermediae-age (e.g., AGB ones) stars in nearby dirr's have shown that some dirr's (e.g., 
NGC 6822 and the Small Magellanic Cloud; SMC) have apparently spherical distributions in 
the outer underlying older stellar populations (e.g., Battinelli et al. 2006; Cioni ct al. 2006; 
Dcmers et al. 2006). These observations combined with the above-mentioned imaging ones 
(e.g., Hunter ct al. 2006) imply that some dirr's possibly have inner main disky structures 
with younger stellar populations surrounded by extended spherical structures composed of 
older stellar populations: some dirr's have two distinct (or "two-component" ) structures. It 
is, however, unclear how such two-component structures were formed in some dirr's (e.g., 
NGC 6822). 

The purpose of this Letter is thus the following two. The first purpose is to introduce a 
classification scheme that divides structures of outer underlying older stellar populations in 
dirr's into the following two: "early-type" (spherical) and "late-type" (disky). The possible 
candidates of dirr's with early- and late-types outer structures are summarized in Table 1 
and the schematic diagram for the proposed structure classification scheme is shown in Fig. 
1. Outer early-type structures in NGC 6822, IC 10, and IC 1613, are confirmed by Battinelli 
et al. (2006), Demers et al. (2004), and BattineUi et al. (2007), respectively 

Tikhonov (2006) found that young stars in a dIrr form a thin disk with the size corre- 
sponding to its visual size whereas older ones form a extended thick disk for a number of 
dirr's (e.g., DDO 216, ESO 381-018, IC 1574, and IC 3104): these dirr's with outer late-type 
structures may well be regarded as normal dirr's. Furthermore Demers et al. (2003) found 
that IC 1613 does not have an outer spherical halo composed of old stars (e.g., red giants) 
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but shows a disky structure composed of carbon stars: IC 1613 is an example of dirr's with 
late- type structures. 

The second purpose of this paper is to demonstrate that gas-rich dwarf-dwarf merging 
can transform dirr's with outer late-type structure into those with early-type ones. We 
mainly describe physical properties of the remnants of dwarf-dwarf mergers with different 

model parameters in order to discuss the origin of dirr's with outer early-type structures. 
It should be stressed here that dirr's with early-type structures discussed in the present 
paper are different from "transition type dwarfs" (e.g., Phoenix and DDO 210) which are 
intermediate-class of objects between dirr's and dEs/dSphs and have HI gas mass of at most 
a few lO^M© (e.g., Grebel et al. 2003): formation and evolution of the transition- type dwarfs 
are discussed in detail by Skillman et al. (2003). We do not intend to discuss dirr's in the 
Local Group, though observations of them provided some clues to formation and evolution 
of dirr's in general (e.g., Mateo 1998). 

Recent observational studies on stellar kinematics in dirr's have shown that there can 
be kinematical differences between stars and HI gas: Hunter et al. (2002) found that stellar 
components in NGC 4449 have a much smaller amount of rotation in comparison with gaseous 
ones. Since we consider that dirr's with these gas-star kinematical differences are possible 
candidates of those with outer early-type structures, we list the candidates in Table 1: such 
star-gas kinematical differences in NGC 6822 and the SMC are observationally confirmed by 
Demers et al. (2006) and Harris & Zaristky (2006), respectively. 

2. The model 

We investigate chemodynamical evolution of gas-rich mergers between dirr's with stellar 

and gaseous disks embedded in massive dark matter halos by using our original Nbody/SPH 
(TREESPH) codes with models of chemical evolution and globular cluster formation (Bekki 
et al. 2002). Since the details of numerical techniques and methods (e.g., ways to model 
chemical enrichment processes) are already given in above paper and our forthcoming papers 
(Bekki & Chiba 2008), we briefly describe them in the present paper. We adopt the Burkert 
profile (Burkert 1995) for the radial density profile of the dark matter halo of a dirr. The 
total masses of the dark matter, the stellar component, and the gaseous one for the dIrr are 
set to be 8.0 x IO'^Mq, 4.0 x lO^M©, and 8.0 x IO^Mq, respectively. 

The dIrr is described as a disk and the radial (R) and vertical (Z) density profile of the 
initially thin are assumed to be proportional to exp{—R/Ro) with scale length Rq — 1 kpc 
and to sech^{Z/Zo) with scale length Zq — 0.2i?o, respectively. The HI diameters of gas-rich 
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galaxies are generally observed to be significantly larger than their optical disks (Broeils & 
van Woerden 1994). Therefore the ratio (sg) of the stellar disk size (rg) to the gaseous one 
(fg) is assumed to be a free parameter. We however mainly show the results of the model 
with Sg=4.0. Star formation in gas is modeled by converting the collisional gas particles 
into collisionlcss new stellar particles according to the Schmidt law (Schmidt 1959) with 
exponent 7 equal to 1.5 (e.g., Kennicutt 1998). The stars formed from gas are called "new 
stars" whereas stars initially within a disk are called "old stars" throughout this paper. The 
simulations have mass and size resolutions of lO^M© and 50 pc, respectively, for stars in all 
models. 

The mass ratio of the two merging dlrr's {1712), the pericenter distance (rp), and the 
eccentricity (cp) are assumed to be free parameters. The orbit of the two dlrr's is set to be 
the same as the xy plane and the distance between the center of mass of the two dlrr's is 20 
kpc. The spin of each galaxy in a merger is specified by two angles 9i and 0^, where suffix i is 
used to identify each galaxy. 9i is the angle between the z axis and the vector of the angular 
momentum of a disk. 0j is the azimuthal angle measured from the x axis to the projection 
of the angular momentum vector of a disk onto the xy plane. Although we run many models 
with different orbits, we mainly show models with Tp = 1 kpc, Cp = 1.0, 9i = 30°, 62 = 120°, 
01 = 90°, 02 = 30°. 

We mainly describe the results of the four representative models with m2= 0.1, 0.3, 
0.5, and 1.0, because the merger remnants of these models show (i) early-type morphologies 
in their underlying older stellar populations and (ii) diversity in the morphologies. We do 
not intend to discuss star formation histories, chemical evolution, and formation of globular 
clusters in dlrr's with outer early-type structures and their dependences on model parame- 
ters: these will be discussed in our forthcoming papers. The stellar-to-mass-to-light-ratio in 
the S-band (Ms/Lb) is assumed to be 1.78 for old stars and 0.25 for new ones so that we 
can convert the simulated stellar mass densities into the -B— band surface brightness (//r)- 
The adopted M^/L-q for old (new) stars are from the tables of stellar population synthesis 
models with [Fe/H] — —1.28 and ages of 5 (0.5) Gyr by Vazdekis et al. (1996). 

3. Results 

Fig. 2 shows that spatial distributions of old and new star are significantly different in 
the sense that spatial distributions of new stars are much more compact and flattened than 
those of old ones in the four representative models with different 1712- Spatial distributions 
of underlying old stars depend on 1712 such that they are more spherical in the models with 
larger m2. Extended gas disks surrounding new stars can be clearly seen in all of the four 
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models and very low- level star formation with the rate of 0.03 — O.IM© yr~^ is still ongoing 
within the gas disks in a sporadic manner for the models. The formation of extended gas 
disks in merger remnants is due essentially to the presence of extended gas disks in merger 
precursor dirr's. Thus these remnants with outer (flattened) spheroids composed of old stars, 
disky structures of young stars, and extended gas disks can be morphologically classified as 
dirr's with outer early-type structures. 

About 30% of initial gas masses corresponding to ~ (2 — 5) x IO^Mq can be still within 
the remnants, which means that the remnants can be regarded as gas-rich dwarfs. Fig. 3 
shows that the gas mass fractions (/g) are significantly higher in the outer parts of the merger 
remnant for the model with m2 = 0.3 (/g ~ 0.6 for ~ 5 kpc). This radial dependence of 
gas mass fraction can be seen in models with different m2, which implies that dirr's with 
early-type structures formed from merging of dirr's have most of their HI gas in their outer 
parts. Fig. 3 also shows that although there is a strong concentration of new stars in the 
central region (i? ~ 1) kpc of the merger remnant (/ns > 0.5), the dark matter halo already 
dominates in mass there (i.e., /b < 0.3). 

Fig. 4 shows that S— band surface brightness distributions (//b) estimated separately for 
old and new stars in the merger remnant with m2 = 0.3. New stars have significantly higher 
than old ones within the central 5 kpc of the remnant so that the optical morphology 
of the remnant can be determined largely by the distribution of new stars: in spite of its 
regular (fiattened) spherical morphology in the underlying old stars, this remnant can be still 
classified as a dlrr rather than a dSph (or dE) in the canonical morphological classification 
scheme, because the outer spherical distribution of old stars {R > 1 kpc) with very low iib 
(> 27 mag arcsec"^) is hard to be detected. It is confirmed that low surface brightness outer 
envelopes composed mostly of old stars are common in the remnants of dwarf-dwarf mergers 
with different m2 in the present study. Thus the present study demonstrates that two dirr's 
with late- type structures can be transformed into one dlrr with an outer early-type structure 
through gas-rich merging. 

4. Discussion and conclusions 

The present study has shown that if gas-rich dirr's have extended gas disks, merging 
between the two can hardly transform them into one gas- free dSph's owing to the presence of 
the gas disk and star- forming regions in the remnant. The Faint Irregular Galaxies GMRT 
Survey {FIGGS) has recently revealed that the medial ratio of rg/rg is about 2.4 for dirr's 
with a median Mb ~ —13 mag (Begum et al. 2008). This observation combined with the 
present results implies that morphological transformation from two gas-rich dirr's into gas- 
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free dSph's via merging is highly unhkely, at least, in the present universe. We thus suggest 
that the remnants of merging between dirr's can be identified as dirr's with outer early-type 
structures and gas such as NGC 6822 (Demers et al. 2006). Possibly, stripping of HI gas 
due to tidal interaction with luminous galaxies and ram pressure of hot gas in group and 
cluster environments could drive further evolution from dirr's with early-type structures into 
gas-free dSph's. 

Given that merging of low-mass galactic building blocks can be a fundamental mode 
of galaxy formation in the hierarchical clustering scenarios (e.g.. White & Rees 1978), the 
present dirr's with outer early-type structures could be the relic of such hierarchical galaxy 
formation. However, it is observationally unclear what fraction of dirr's are actually those 
with early-type structures owing to the lack of extensive statistical studies on spatial dis- 
tributions of older stellar populations in the outer parts of dirr's. van den Bergh (1988) 
suggested that typical dirr's are triaxial systems with axis ratio of 1.0:0.9:0.4 based on the 
observed distributions of inclination angles of nearby dirr's. Observational studies on the 
distribution of ellipticities in light distributions of 30 low-surface brightness dirr's suggested 
that the intrinsic shape of the dirr's is triaxial and slightly less spherical than dEs (Sung 
et al. 1998). Although these observations imply that inner main components of dirr's have 
more flattened distributions, it remains unclear whether outer older ones also have such 
distributions. 

As shown in Fig. 4, the outer Tindcrlying structure composed of older stars in a dirr 
formed by merging is hard to be detected owing to its very low surface brightness (e.g., 
liB > 30 mag arcsec"^ for i? > 3 kpc), even if the mass fraction of the older stars is significant. 
This suggests that observations based on spatial distributions of individual bright stars (e.g., 
RGB/AGB/carbon stars) rather than on broadband imaging would be better to reveal the 
presence of outer extended structures composed mostly of older stars in dirr's. Although 
the presence of intermediate-age stars within the outer halos for a number of dirr's has been 
reported (e.g., Albert et al. 2000), the details of their spatial distributions are not so clear 
for some of the dirr's. As suggested in the present simulations, the possible presence of two 
distinct structures (i.e., inner disks and outer spheroids) would have valuable information on 
the formation of dirr's. Thus it is doubtlessly worth while for future observations to detect 
outer spherical or disky structures composed of older stars in dirr's. 

Although nearby dirr's are observed to show signs of rotational kinematics in their 
gaseous components (e.g., Mateo 1998 for a review), it is observationally unclear whether 
their stellar components have also rotational kinematics (e.g.. Hunter et al. 2002; 2005). 
The present study predicts that V/a can be significantly different between gas and stars 
in the merger remnants in the sense that their stellar V/a are smaller than gaseous ones. 
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Therefore, gas-rich dirr's with kinematically hot stellar components are promising candidates 
of the remnants of mergers between dirr's with extended gas disks. Thus future systematical 
studies of stellar kinematics of dirr's will help us to understand what fraction of dirr's are 
those with early-type structures and thus possibly formed from dwarf-dwarf merging. 

Both the presence of intermediate-age stars in very outer regions of dirr's (e.g., Albert 
et al. 2000; Letarte et al. 2002) and apparently spherical distributions of these stars (e.g., 
Battinelli et al. 2006) can be explained by the present merger scenario of dirr formation, if 
merging can happen 2 — 10 Gyr ago. However, dirr-dirr merging can be only one of several 
possible physical mechanisms for the formation of outer structures composed of intermediate- 
age stars. For example, minor merging of very tiny dwarfs with intermediate-age stars onto 
a dIrr can be also responsible for the presence of AGE and carbon stars in the outer region 
of the dIrr, though such minor merging can hardly produce an outer spherical structure. If 
supersonic gaseous outflow driven by central star-formation activity in a dIrr can interact 
with the outer halo gas and consequently trigger star formation in the high-density shocked 
gaseous regions, such new stars might well be identified as intermediate-age halo stars several 
Gyrs later: dynamical relaxation after the formation of halo stars can be responsible for the 
formation of a spherical structure around the dIrr in this scenario. We thus suggest that 
physical properties of older and intermediate-age stars in the very outer regions of dirr's can 
provide valuable information on formation and evolution of dirr's. 

I (K.B.) am grateful to the anonymous referee for valuable comments, which contribute 
to improve the present paper. K.B. acknowledges the financial support of the Australian Re- 
search Council throughout the course of this work. The numerical simulations reported here 
were carried out on SGI supercomputers at Australian Partnership for Advanced Computing 
(APAC) in Australia for our research project. 
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Table 1. Possible candidates of dirr's with outer early-type (i.e., spherical) structures, 
outer late-type (i.e., disky) ones, and kinematical differences between stars and HI gas. 



dIrr properties 



Examples 



Outer spherical structures 
Outer disky structures 
Star-gas kinematic differences 



NGC 6822, IC 10, IC 1613 

DDO 216, ESO 381-018, IC 1574 
NGC 6822, SMC, NGC 4449 



- 10 - 



Face-on Edge-on 



Late-type 



dirr 



Early-type 





# Old stars 

# New stars 
O HI gas 




HI stripping 

dSph/dE ? 



Fig. 1. — A schematic diagram for a possible diversity in outer underlying structures of 
older stellar populations (e.g., RGB/AGB stars) in dirr's. Here dirr's are divided into (i) 
those with outer "early-type" (i.e., spherical) structures and (ii) those with "late-type"' (i.e., 
disky) ones according to the outer underlying structures of older stellar populations. Old 
stars, new (or young) ones, and HI gas are shown in red, blue, and yellow, respectively, both 
for the face-on and the edge-on views, dirr's with late-type structures (e.g., DDO 216) are 
basically dwarfs with rotating stellar and gaseous disks whereas those with outer early-type 
structures (e.g., NGC 6822) have both spherical structures of older stars and disky ones 
of young and (old) stars and HI gas. The degree of flattening can be different for outer 
spheroids and the spheroids can be supported dynamically by velocity dispersion and by 
rotation. Possible stellar halos composed of very old stars (e.g. Minniti & Zijlstra 1996; 
Aparicio & Tikhonov 2000) are not used for this morphological classification. It should 
be here stressed that owing to the presence of intermediate-age stars, the outer spherical 
structures are different from stellar halos dominated by very old stars with ages larger than 
10 Gyrs in luminous galaxies like the Galaxy, dirr's with outer early-type structures are 
different from "transition type" (dSph/dIrr) dwarfs owing to the presence of star- forming 
regions and abundant HI gas. One of possible evolutionary links between dirr's with late- 
type and early-type structures is shown: dirr's with late-type structures can evolve into blue 
compact dwarfs (BCDs) and further into dirr's with early-type ones by dwarf-dwarf merging 
(Bekki 2008). Complete or partial stripping of HI gas via some physical processes (e.g., 
tidal interaction with other galaxies or ram pressure) would transform dirr's with early-type 
structures into dIrr/dSph transition dwarfs and gas-free dSph's (or dE's). 
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Fig. 2. — Final mass distributions of old stars (magenta), gas (cyan), and new stars (yellow) 
in the remnants of dwarf-dwarf mergers for the four representative models with m2=0.1 (a), 
0.3 (b), 0.5 (c), and 1.0 (d). A bar in the lower left corner in each panel measures 2 kpc. 
It is clear that the remnants have inner disky (or very flattened spherical) structures of new 
stars, outer (flattened) spherical ones of old stars, and extended HI gas disks surrounding 
new stars. These remnants can be classified dirr's with outer early-type structures. 
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Fig. 3. — The radial dependences of mass fractions in the merger remnant for the model 
with 7722 = 0.3: /b is the mass fraction of baryon (gas and stars) to all components (solid), 
/g is that of gas to baryonic components (dotted), and f^s is that of new stars to all stellar 
components (dashed). 




Fig. 4. — The 5-band surface brightness profiles (fi^) for old stars (thin) and new ones 
(thick) for the merger remnant in the model with m2 = 0.3. 
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